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SUMMARY 


A vector study of the partial-differential equation of steady 
linearized supersonic flow is presented. General expressions, which 
relate the velocity potential in the stream to the conditions on the 
disturbing surfaces, are derived. In connection with these general 
expressions the concept of the finite part of an integral is discussed. 

A discussion of problems dealing with planar bodies is given and 
the conditions for the solution to be unique are investigated. 

Problems concerning nonplanar systems are investigated, and methods 
are derived for the solution of some simple nonplanar bodies. The sur- 
face pressure distribution and the damping in roll are found for rolling 
tails consisting of four, six, and eight rectangular fins for the Mach 
number range where the region of interference between adjacent fins 
does not affect the fin tips. 


INTRODUCTION 


In the presentation of the theory of the flow of an idealized 
incompressible fluid, vector methods can be used to reduce greatly the 
mathematical manipulations involved. The study of steady linearized 
supersonic flow may also be aided by the use of vector methods. Two 
types of approaches, however, can be used. Perhaps the more obvious is 
to make use of common vector methods as was done in reference 1. The 
other vector method, which was introduced by Robinson in reference 2 
and is used in this paper, appears to be more suited to the study of 
the linearized partial-differential equation of steady supersonic flow. 
This method allows a derivation of a hyperbolic scalar potential and a 
hyperbolic vector potential along lines analogous to the derivation 
sometimes used (reference 3* ch. VIII) in dealing with common scalar 
and vector potentials. 
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The present paper presents a vector derivation of many general 
results which have been found by various methods and are given in the 
published literature on the linearized partial-differential equation 
of supersonic flow and also presents some results which are not found 
in the literature . The general results of Hadamard (reference 4, 
p. 207), Puckett (reference 5), and Heaslet and Lomax (reference 6) 
are found as special cases of a general expression for a scalar poten- 
tial, and the results found by Robinson (reference 2) are obtained by 
the use of a vector potential. The derivation of the scalar potential 
doubtlessly helps to clarify the concept of the finite part of an 


A discussion of problems dealing with planar bodies immersed in a 
supersonic flow is given, and the conditions necessary for the solution 
to be unique are investigated. 

Problems dealing with nonplanar systems are also discussed, and 
methods are derived for the solution of some simple problems dealing 
with nonplanar bodies. The surface pressure distribution, the span- 
wise loading, and the damping in roll are found for rolling tails con- 
sisting of four, six, and eight rectangular fins for the Mach number 
range where the region of interference between adjacent fins does not 
effect the fin tips. 


integral . 


SYMBOLS 


A 


hyperbolic vector potential 


A 


aspect ratio of tail fin 


a 


positive constant 



arbitrary constants 


c 


chord 


f 


arbitrary vector functions 

scalar function defined by equation (19) 


G 


vector function associated with vector function F 


H 


vorticity vector 


b/2 


span of tail fin 
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i, 

M 

n 

n h = -i(B 2 V 1 + 

v - -«>V 
v - - 1 |j 2 v i* 


unit vectors in x-, y-, and z-directions, respectively 
Mach number 

unit vector normal to element of area da 

jv 2 + 

+ jv 2 * + kv 3 ’ 

. * , * 

+ JV 2 + kv 3 


P 

Q 

q 

q’ 

r = 

R' 


pressure-difference coefficient 
rate of roll 

function used in equation of surface of discontinuity 

part of velocity vector which is made up of hyperbolic 
curl of vector potential 

total perturbation velocity 

/(x - |) 2 - 3 2 (y - ti) 2 - p 2 (z - £) 2 
small constant 


r = ]j(x - 4) 2 
S 

S 0 

s 5 ,s 6 ,t,t- 

V 

v 0’ v ±> v 2 


+ pNy - t ) 2 + p\ z - o 2 

area of tail fin 
surface of discontinuity 

surfaces of integration 

free -stream velocity 
volumes of integration 


W = | £ 


x,y,z 


Cartesian coordinates (x-axis parallel to free- stream 
direction) 
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= n/m 2 - 


spanwise circulation 
small positive quantity 


' r T.E. .. p T.E. 

2 f 0 X dx = | r 
^ Jl.e. ^ Jl.e. 


ACp dx 


e,p 

X 


v l ,v 2 ,v 3 
V 1 > v 2 > v 3 


v l*> v 2*' v 3* 


Cartesian coordinates (£-axis parallel to free-stream 
direction) 

polar coordinates 

scalar functions 

given volume 

direction cosines of outward normal to element of area da 

direction cosines of normal (directed away from point 
(x,y, z) ) to surface S Q 

direction cosines of normal to element of area da used 
in equation (45) 


a slope of deflected area 

t area of integration 

scalar Potentials 


rolling-moment coefficient per fin 

( Rolling moment per fin\ 

I pv2 s| / 


c z = 

6 P 


dC- 




p — > 0 


indicates integration over closed line or surface 


denotes finite part of integral 
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THEORY 


This paper deals with the linearized partial-differential equation 
of steady supersonic flow. This equation is given by 

-(3 2 0 + 0 + d> = o (l) 

xx r yy *zz (L) 

The potential is assumed to be continuous in the stream direction, and 
the potential is assumed to be always finite. Assuming the potential 
to be finite and continuous in the stream direction has the effect of 
requiring the aerodynamic lift and moment (calculated by use of the 
linearized pressure) of finite bodies to be finite since the linearized 
pressure is related to the derivative of the potential in the stream 
direction. The expression linearized pressure" refers to the pressure 
obtained by neglecting all powers of the perturbation-velocity compo- 
nents above the first. 


Vector Operators and Identities 

Certain operators, which are closely associated with the linearized 
hyperbolic partial- differential equation of supersonic flow (the two- 
dimensional wave equation) , are added to the vector operators commonly 
used. The basic operators have been used previously in references 2 
and 7- 

The gradient operator is defined by 

. 6 .6 S 

V=1 Sx + J Sy + k ^z 

The analogous hyperbolic gradient operator defined by Robinson in 
reference 2 may be expressed as 


vh _ -i|3 2 JL + j JL + k 

The hyperbolic divergence of an arbitrary vector E is given by 

Vh • E 

Similarly, the hyperbolic curl of the vector E is given by 


Vh x E 
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The divergence of the gradient operator is sometimes denoted by 

2 a 2 a 2 a 2 

V = V * V = — p + — p + — p 
6x dy 6z 

The analogous divergence of the hyperbolic gradient operator is denoted 

a. = - . ^ a 2 . a 2 


V h 5 V • Vh =-3 


dx 2 3y 2 6z 2 


The following identities are needed. Let E be a vector and and A 
be scalar functions of x, y, and z. Then, 


Vty • VA = VA • 

(2a) 

V • tyE = tyV ' E + E • Vty 

(2b) 

V X (V X E) = V(V * E) - V 2 E 

(2c) 

V • (V X E) = 0 

(2d) 

Vhi • VA = VhA • Vt|r 

(2e) 

Vh • tyE = \^Vh • E + E • Vh^ 

(2f) 

V x (vh X E) = Vh(v • E) - V 2 hE 

(2g) 

Vh X (V X E) = V (Vh • E) - V 2 hE 

(2h) 

Vh • (vh X E) = 0 

(21) 


These identities can be proved by direct expansion. 
The divergence theorem may be expressed as 


E • n da = V • E dv 


( 3 ) 


where n is the normal unit vector to the element of area da. The 
vector n is expressed mathematically as 

n = iv 1 + jv 2 + 
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where and are the direction cosines of the outward drawn 

normal to the element of area da. 


A theorem more general than the divergence theorem is given hy 
(this theorem follows from the results of reference 8 , p. 87 ) 


§ (w 


+ CgVgEy 


+ C 3 v 3 E z J 


da = 



c*E v dE z \ 

+ c 2 + c 3 3T-) dv 


where_the subscripts x, y, and z refer to components of the vec- 
tor E, and C-^ Cq, and C 3 are arbitrary constants. Note that if 

= C 2 = C 3 = 1 the preceding equation reduces to equation (3). If 

c 1 = -P 2 



= 1 


the preceding equation reduces to 




* v 2 E y + v 3 \) 


da 



Sr + 5F + 



dv 


or 


/ E • 


da = / Vh • E dv 


w 


where 


n h = -i P^ v i + ^ v 2 + ^3 

If the divergence theorem as expressed by equation (3) is applied 
to a volume throughout which 

V • E = 0 

then the surface integral over the bounding surface is 
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provided that no surfaces exist_inside the volume of integration across 
which the normal component of E is discontinuous. Similarly, if 
equation (4) is applied to a volume throughout which 


Vh • E = 0 


then the surface integral over the hounding surface is 



n^ da = 0 


provided that there are no surfaces inside the volume of integration 
across which E • is discontinuous. It is interesting to note, 

however, that surfaces exist inside the volume of integration across 
which E • n can he discontinuous while at the same time E • n^ 
remains continuous. It follows that for such a surface n and 
must satisfy the relation 


n • n h = 0 (5) 

Let Q(x,y,z) =0 he the equation of such a surface. 

Then , 


and 


n 



VQ 


n- 


h 



where the subscripts indicate differentation. Substituting the 
preceding expressions for h and n^ into equation (5) yields 


-3 


2 



= 0 


( 6 ) 


Any solution of equation (6) set equal to zero is the equation_of a 
surface across which V * E may be discontinuous while Vh • E remains 
continuous. The fact that the Mach cone from any arbitrary point satis- 
fies equation (6) can be easily verified. The equation of the envelope 
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of the Mach cones from an arbitrary line also satisfies equation (6) 
(reference 9 , p. 106). 


In the following sections use is made of the concept of the finite 
part of an infinite integral. This concept was introduced by Hadamard 
(reference 4) and has been used by a number of other investigators. 

The finite part is, however, sometimes confusing. This section was 
therefore included in an attempt to give a realistic picture of the 
finite-part concept .and also to present the first steps of the deriva- 
tion of the scalar and vector potentials. 

The concept of the finite part of double integrals as defined by 
Hadamard and used in this paper is different from the concept of the 
finite part of double integrals as defined in reference 10. The essen- 
tial difference between these two definitions lies in the manner in 
which the singular points along the Mach cone are treated. 

In reference 3, page 183, a vector function is used in the deriva- 
tion of the common scalar and vector potentials. The analogous vector 
function based on equation (1) is 


The hyperbolic divergence of vector W with respect to variables 
I, T), and £ is given by 


The preceding equation indicates that the hyperbolic divergence of the 
vector W set equal to zero yields the partial-differential equation 
°f linearized supersonic flow. A mathematical derivation of W can be 
obtained; however, for the purposes of this paper such a derivation is 
not needed. 


Finite Part of Integrals Which Arise 


in Steady Supersonic Flow 


w = gV0 -Ifci 


where 


R = - ^) 2 - P 2 (y - n ) 2 - 3 2 (z - £) 2 


Vh • W 
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The result of applying equation (4) to the vector W is 


$ * ' 4 " g) 


. da = 


I 


R 


i V^h0 dv 


( 7 ) 


When 0 satisfies equation (l) throughout the volume of integration, 
the right side of equation (7) is zero) thus. 


Equation (7) is applied to a volume (denoted by v Q ) enclosed in 
the forward Mach cone from the point (x,y,z). This volume is bounded 
by the surface given by R = R 1 , where R' is a small constant, and an 
arbitrary surface Sq enclosed in the forward Mach cone from the point 
(x,y,z). A cross section of the region of integration is shown in fig- 
ure 1. Rote that this region is analogous to the region that is some- 
times used in calculating the potential function satisfying Laplace's 
equation (reference 3, pp. 151-153)- For regions such as the one shown 
in figure 1, equation (7) may be written as 



( 8 ) 


when 


V 2 h <f> = 0 



. n^ da + 



. n^ da = 



(9) 


where T represents the area of integration when R = R' . 


The integral over the area T may be reduced to 



( 10 ) 


where r is given by 



(ID 
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Since R' is a constant, equation (10) can be written as 



Equation (9) can now be written as 

( v <* • 5 h * ^ (s ^ ■ "h 38 =/ iv 

( 12 ) 

If 0 is required to satisfy the linearized partial-differential 
equation of steady supersonic flow, then 

V 2 h0 = 0 

and equation ( 12 ) reduces to 

( V ^ ' % + ^) da + J ^ |) • % da = 0 ( 13 ) 

S 1 ' 

If R’ is made smaller and smaller the integrand of the integral over 
the area T in equation ( 13 ) remains finite except on the small area 
close to the point (x,y,z). In anticipation of taking the limit of 
equation (13) as R' approaches zero, the small area close to the 
point (x,y,z) is removed from the area T. The area T is divided 
into two parts. One part is the area of T which is downstream of the 
surface given by 


£ = x - € 

where e is small but larger than R' . This area is denoted by t. 
The remaining part of T (denoted by T') is the area of T which is 
upstream of the surface 


I = x - e 

A cross section of the region of integration with T divided into t 
and T' is shown in figure 2 . Equation (13) can now be expressed as 
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R' 






da + 


J /i - 07 . n h da = 0 

S 1 N 

where R T is smaller than e. 


(14) 


Since 0 is continuous and therefore its values over t are 
approximately constant for small values of e , the integral over the 
area t can be written as 



n^ da + 


(3 2 0(x,y,z) 

R' 



da 

r 


(15) 


When the second integral of expression (15) is integrated, equation (14) 
becomes 



n h 


2*e0(x,y,z) 

da + 





da = 0 


(16) 


If R' is made to approach zero, equation (l6) applies even to the 
limit where R f is zero. 

The limit of equation (l6) as R’ approaches zero may he written 
as 


lim 


R' 


k 0 



2ng0(x,y,z) 

R’ 


2«0(x,y,z) 




(17) 


The integrands of the integrals over the areas t and T T are always 
finite and it can he shown that their first derivatives with respect 
to R* approach zero as R’ approaches zero; therefore, the product 
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of 1/R' and these integrals either approaches zero in at least the 
order of R' or approaches infinity as R' approaches zero. Thus it 
follows that the integrals over the areas t and T' have no finite 

terms remaining after the limit (R' > 0) has been taken. The sum of 

the terms of equation (17) must be zero; thus the singularities 
resulting from the integrals over the areas t and T' must cancel 
the singularities which arise from the integral over the area Sq. 

From the preceding considerations it follows that one method of 
evaluating the finite part of infinite integrals of the type appearing 
in equation (17) is to evaluate the integral when R' is small but not 
zero and neglect the terms multiplied by powers of l/R'. Other infinite 
integrals sometimes arise, however, for which the finite part cannot be 
obtained by neglecting powers of l/R' . For example, if equation (17) 
is differentiated with respect to one of the variables (x, y, or z) 
an equation containing the velocity component is obtained. In some 
cases, when the point (x,y,z) lies on the surface Sp the infinite 
terms are of the order (in R')/R' and of the orders (l/R') n . In these 
cases, the finite part of the infinite integrals can be obtained by 
evaluating the integrals when R’ is small and neglecting the terms 
multiplied by powers of l/R' and (in R')/R'. 

The process of removing the infinite parts of an integral, however, 
has been derived by Hadamard (reference 4, book III, ch. I). Hadamard 
used his methods of evaluating the finite part of integrals in finding 
solutions to certain hyperbolic equations including the linearized equa- 
tion of steady supersonic flow. Perhaps a fact worth noting is that 
the integrals of equation (17) are double integrals and when the methods 
given by Hadamard are used the methods given for multiple integrals 
should be used. In the past, the singular points (points on the Mach 
cone where the derivative of (x - |) 2 _ p 2 (y _ t)) 2 - 3 2 (z - £) 2 with 
respect to the variable of integration is zero) have caused some con- 
fusion; as Hadamard points out (reference 4, p. 147), these singular- 
points must be removed from the area of integration before the finite 
part is taken. Particular attention should be given to paragraph 92 
of reference 4 since the special type of integrals discussed therein 
sometimes arises in dealing with planar problems. 

Robinson (reference 2) has shown that when using Hadamard's methods 
the order of integration may be changed without affecting the finite 
part and that it is permissible to differentiate under the integral 
sign of a multiple integral without considering the variable limits 
which lie along the boundary where the integrand is singular, provided 
that only the finite part is taken. Both Hadamard and Robinson have 
shown that in differentiating an improper integral which has an integrand 
that has a one -half power singularity along variable limits the variable 
limits may be neglected provided the finite part of the resulting 
integral is taken. 
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The term "finite part" is somewhat misleading since the finite 
part of an integral can he infinite. In certain cases the integral 
is infinite even after the terms which approach infinity as R' 
approaches zero have been neglected. 


Scalar Potential 

The preceding arguments show that the finite parts of equation (17) 
can be equated to zero; thus, 


-2«0(x,y,z) 



n^ da = 0 


where the symbol f before the integral denotes that only the finite 
part is to be taken. The preceding equation may be solved for the value 
of the potential at the point (x,y,z); the result of this operation is 
given by 

0(x,y,z) = ^ j (| V0 - p f) . n h da (l8) 

S 1 

It should be remembered that surfaces can exist inside the forward Mach 
cone from the point (x.y.z) across which i V0 - 0V i can be discon- 

- V0 - 0V g- 1 . h h remains continuous . 

Equation (l8) is an expression for the scalar potential at the 
point (x,y,z) in terms of the potential and its derivatives with respect 
to n^ on the surface Sp. A more general expression for the scalar 
potential than that given by equation (l8) can be obtained. If within 
the volume enclosed by the forward Mach cone from the point (x,y,z) and 
the surface Sp equation (l) is not satisfied and V 2 h operating 
upon 0 yields 


= f(s,n,£) 

then equation (10) becomes 


(19) 




n h da 



R 


dv 


( 20 ) 
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Provided that f(t,T|*£) is always finite, the right side of equation (20) 
is finite; furthermore, the right side of equation (20) remains finite 
as R' approaches zero. If in equation (20) R' is made to approach 
zero and only the finite parts of the integrals are retained, then the 
resulting expression is 

?S(x,y,z)=i f (iv0 - 0V i) • TT h da - ±J hMlii dT ( 21 ) 

S 1 V 1 

where vq represents the volume Vq when R' is equal to zero. Equa- 
tion (21) is equation ( 58 ) of reference 4 where p 2 has been set equal 
to one. Note that the volume integral in equation (21) has the appear- 
ance of the integral for the potential resulting from a volume distribu- 
tion of sources in an incompressible flow. 

The assumption has been made that 0 is continuous throughout the 
volume vq. It is also assumed that no surfaces exist inside vq 
across which a0/6n h is discontinuous. If equation (21) is applied to 
a volume v, , which has surfaces across which 0 and/or the derivative 
of 0 in the direction of n h is discontinuous, these surfaces of dis- 
continuity can be removed from the volume of integration by allowing the 
arbitrary surface Sq to envelop them (see fig. 3) • For volumes of 
integration where the surfaces of discontinuity have been removed in 
this manner, the scalar potential can be written as 



where Sq denotes the surface of discontinuity, and A 0 is the poten- 
tial difference across the surface Sq. The notation d/dnjj' is used 
to denote the operator 



a 

sc 


V 


a 

3 ^ 


+ V 


, a 

3 sr 


For the cases where no surface of discontinuity exists inside the 
volume and 0 and V0 are zero on the surface - Sq , equa- 

tion (22) reduces to 
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0(x,y,z) = F dv (23) 

V 1 

From equation (19) 

V 2 h0(x,y,z) = f(x,y,z) (24) 

Note that equation (23) is a solution of the partial -differential 
equation (24). 

For most problems in linearized supersonic flow, f(£,T|,£) is zero 
and f is zero upstream of the disturbing body. For such problems, the 
surface Sq - Sq can be taken to be located upstream of the disturbing 
body where 0 and V0 are zero. In this case, equation (22) reduces 
to 


0(x,y,z) 



1 d j j 6 

R A 0 - A 0 



da 


(25) 


If the surface Sq is confined to the £ = 0 plane, equation (25) 
reduces to equation (10) of reference 6. In this reference the boundary 
conditions for airfoils are discussed. 


Components of Vector Field 

Let F be a vector which is finite and integrable in a given 
volume (denoted by X) and is zero outside the volume X. To each 
point in the volume associate the vector 


G(x,y,z) 



(26) 


where v 2 denotes the part of the volume X enclosed in the forward 
Mach cone from the point (x,y,z). 

From equations (24) and (26), it follows that each component of G 
satisfies the relation 


V 2 hG i (x,y,z) = -2«F i (x,y,z) 


( 27 ) 


where the subscript i refers to any component of the vector field. 
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Let \|r Q (x,y,z) be a scalar and 
the equations 


A(x,y,z) be a vector defined by 


2nt 0 = Vh • G = 


X 



( 28 ) 


and 


f 


X 


2nA = V X G = 


XV^dv 


(29) 


Equation (2h) indicates that 


Vh x (v X G) = V(Vh • G) - V^hG 


(30) 


Substituting the expressions for V x G, Vh • G, and V 2 hG as given 
by equations (_27), (28), and (29), respectively, into equation (30) and 


Since F is an arbitrary vector, equation (31) indicates that any 
finite integrable vector field can be expressed in terms of the gradient 
of a scalar and the hyperbolic curl of a vector. Equation (31) has the 
appearance of the Helmholtz theorem (reference 3, p. l87)j however, 
since Vq and A are found by integration only in the forward Mach 
cone from the point (x,y,z), equation (31) hardly seems to he a state- 
ment of the Helmholtz theorem as is commonly given. The result given 
by equation (31) was obtained by Robinson in reference 2. 


Equation (31) indicates that the perturbation velocity vector can 
be divided into two parts. One part is the gradient of a scalar func- 
tion, and the other is the hyperbolic curl of a vector function. The 
vector function is analogous to the common vector potential (reference 3 
pp. 104 and l88)j therefore, the vector function is referred to as the * 
hyperbolic vector potential. Thus, if q' is the total perturbation 
velocity vector, then 


solving for F yields 


F(x,y,z) = -Vt 0 (x,y,z) + Vh x A(x,y,z) 


(3D 


Hyperbolic Vector Potential 


q' = V0 + Vh X A 


(32) 
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where 0 is the scalar potential and A is the hyperbolic vector 
potential. The part of the velocity vector which is made up of the 
hyperbolic curl of the vector potential is denoted by q. 

By direct expansion it can be shown that 

Vh • q' = Vh • V0 + Vh •• (Vh X A) = V 2 h0 + Vh • (Vh X A) = 0 (33) 

Equation (33) indicates that the hyperbolic divergence of the perturba- 
tion velocity vector is zero. 

The vorticity vector is given by 

H = V X q* (34) 


Therefore, from equation (32), 

H = V X (Vh X A) 


or 

H = Vh(V • A) - V 2 hA 

From equations (2d) and (29), the divergence of the hyperbolic 
vector potential is zero; thus, 

H = -V^hA (35) 

Each component of equation (35) is a partial-differential equation 
of the form of equation (24) j thus, from equation (23) each component 
of equation ( 35 ) has a solution given by 

x i r H i^^) 

A i (x,y,z) = / r dv ( 36 ) 

J n 

where the subscript_ i refers to any component of the vector H. Since 
each component of A is given by equation (36), then 

A(x,y,z) = d r dv (37) 

V 1 

The velocity vector resulting from the hyperbolic vector potential is 
therefore given by 
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or 


q = Vh X A = ^ / Vh x dv 

J V]_ 


(38) 


u = 


5 2 f r (y - Ti)H z - (Z - 5)h 3 


2jt 


dv 


(39a) 


,2 *n (z - g)H x - (x - ?)H. 


V = 


2 * 


dv 


(39b) 


? 2 f C (x - 5)H V - (y - rj)lL 


w = 


2n 


dv 


(39c) 


where the subscripts refer to the components of the vector H. The 
results given by equations (39) were obtained by Robinson in reference 2. 


Vortex Sheets 


If the vorticity is confined to a surface S 2 , equation (37) 
becomes 


A(x,y,z) = i J da (40) 

Equation (40) is an expression for the hyperbolic vector potential 
resulting from a surface of vorticity. Note that if the vorticity is 
zero except on the surface S 2 , then equation (35) reduces to 

V 2 hA = 0 


By removing_the surface S 2 from the volume of integration each com- 
ponent of A can be expressed as (from equation (25)) 


Ai(x,y,z) 


1_ 
2 it 



d V 


- AAi 


6 



(41) 


’2 
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where the subscript i refers to any component of the vector A. 

Since each component of A is given by equation (4l), then 

A(w)"S t / S 2 (E i ^- M '^f) 4 ' (42) 

Note that if AA is zero, equation (42) reduces to 

*U,y,*) - S A SS7 da (43> 

u S 2 h 

By comparing equations (40) and (43) it follows that, on the surface S 2 , 

h = m 

n 

Equation (44) indicates that across a surface of vorticity the 
derivative of the hyperbolic vector potential in the direction of n^' 
is discontinuous. Thus, a lifting surface can be represented by a 
continuous hyperbolic vector potential, while it can be shown that a 
thickness effect can be represented by a discontinuous hyperbolic 
vector potential. Note the contrast with the scalar potential, which 
uses a continuous potential to represent a thickness effect and a dis- 
continuous potential to represent a lifting surface. 


Further Development of Scalar Potential 


The scalar potential can be expressed in forms other than those 
already presented. Equation ( 8 ) is applied to the region bounded by 
the arbitrary surface Sp, the forward Mach cone from the point (x,y,z), 
and a second arbitrary surface S 3 enclosed in the forward Mach cone 
from the point (x,y,z) and upstream of the surface- Sp. A cross section 
of such a region is shown in figure 4. The result of applying equa- 
tion ( 8 ) to this region is 



1 80 ’ 

R 3nJ* 


- 0 ' 


<3nw* 



60' 


- 0’ 




||da 


0 


(45) 


provided that 0' is a solution of equation (l). Note that the scalar 
potential as given by equation ( 18 ) is independent of 0 ' so that 0 ' 
is arbitrary so long as it satisfies equation (l) throughout the proper 
volume . 
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If for a finite distance upstream 0' is zero and remains zero 
for greater distances upstream, the surface S 3 may he chosen in this 
region so that the integral over S 3 in equation ( 45 ) is zero; thus. 


2 * 



* i) 

5n^* R J 


da = 0 


(46) 


Equations (l 8 ) and (46) can be combined to yield 


0 (x,y,z) 




- (0 - 0 ’) 


6 1 
8 n h * R 


da 


(47) 


The only restrictions placed on 0 r at this point are that it satisfy 
equation (l) and be zero at a finite distance upstream. In many cases 

0 may be chosen so that ^ % or 0 - 0 T is zeroj therefore, 

in these cases, 0 can be expressed as 


0(w) 


da 


(48a) 


or 


0 (x,y,z) 



^ - *' >3^ S ds 


(48b) 


Equations (48) are quite useful; however, remember that they apply only 
when 0 ’ can be chosen so that 0 ! does not violate any of its 
restrictions . 

equations (48) can be applied to problems where either 0 
or 00 /dnh is given on the surface Sq. The application of these 
equations to most nonplanar problems of either type, however, lead to 
quite unwieldy integral equations. 
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APPLICATIONS 
Planar Problems 


Many problems in linearized supersonic flow deal with the surface 
of discontinuity confined to a plane surface parallel to the x-axis. 

In this section a general discussion of this type of problem is given. 
The coordinates are located so that the surface of discontinuity is in 
the £ = 0 plane. 

The scalar potential at an arbitrary point (x,y,z) above the 
£ = 0 plane is (from equation (47)) 


0(x,y,z) 



(h-*n)s+V 



d! dq 


(**9) 


In this case, the surface Sp is the £ = 0 plane. 

If 0'(x,y,O) is chosen equal to 0(x,y,O) the potential becomes 
symmetric with respect to the £ = 0 plane. Thus, for z = a 


0(x,y,a) = 0' (x,y, -a) 


and 


0 z (x,y,a) = -0 z '(x,y,-a) 
For this case, equation (^9) reduces to 


0(x,y,z) 


l f 0 Z (6,T),O) 
rt J R 

b l 


d£ dn 


(50) 


Equation (50) was given by Puckett in reference 5» 

If 0'(x,y,O) is chosen equal to -0(x,y,O) the potential 
becomes antisymmetric with respect to the £ = 0 plane. In this case, 

0(x,y,a) = -0' (x,y, -a) 


0 z (x,y,a) = 0 Z ' (x,y, -a) 


and 
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Thus, equation (49) is reduced to 


0 (x,y,z) 



R 3 


d£ dq 


(5D 


Note that for surface Sp not confined to a plane parallel to the 
x-axis, a choice of 0'(x,y,z) at the surface Sp to equal 0(x,y,z) 

^ 2 _(x>y>z) a p sur f ace Sp to equal 

d0' (x v z) dn h 

’ ’ at Sp. Similarly, choosing 0'(x,y,z) at the surface 

.al to -f> 

c>0(x,y , z) 


at Sp does not cause 

cSv 

Sp equal to -0(x,y,z) at Sp does not cause d0' (x,y>z) at Si to 


equal 


3n 


at Sp. 




Provided the discontinuities are restricted to the £ = 0 plane, 
the scalar potential can also he expressed as follows (from equation ( 18 )) 


0(x,y,z) =-i 


5 — + 0(l,n>o)^- 

R J 


dl dq 


(52) 


for positive z. A comparison of equations (50) and ( 51 ) with equa- 
tion ( 52 ) shows that the two terms of the integrand of equation ( 52 ) 
contribute equal amounts to the potential at any point (x,y,z). 


Since the terms of the integrand of equation (52) contribute equal 
amounts to the potential at the point (x,y,z) as z approaches zero, 
equation ( 52 ) must reduce to 


0(x,y,O) = -~ 


0z(£> T b o ) 


R 


d| dq 


0 (x,y,O) 


The preceding equation can also be obtained by examining the limit of 
equation (52) as z approaches zero. If this procedure is done the 
entire contribution of the second term of the integrand of equation ( 52 ) 
is found to come from the point at the apex of the hyperbola formed by 
the intersection of the Mach cone from the point (x,y,z) and the 
£ = 0 plane. Note that if the integration is performed first with 
respect to q then, when the methods of Hadamard are used, the point 
at the apex of the hyperbola is a singular point and must be removed 
from the area of integration by a process such as is given in refer- 
ence 4, page 147. 
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If 0 Z is prescribed over the £ = 0 plane, then the potential 
is given uniquely by equation (50). Similarly, if the potential is 
prescribed over the £ = 0 plane, then the derivative of 0 with 
respect to z is determined over the z = 0 plane. This result fol- 
lows from equation (51) since prescribing 0 over the z = 0 plane 
determines the potentials in the space above the z = 0 plane; there- 
fore, it also determines 0 Z in the space above the z = 0 plane and 
the limit of 0 Z as z approaches zero from the positive direction. 

The question that arises is whether 0(x,y,z) is uniquely deter- 
mined in the space above the z = 0 plane if 0 is prescribed over 
certain areas of the z = 0 plane and 0 Z is prescribed over the 
remaining areas. If the assumption is made that 0 is not determined 
uniquely, then at least two potential functions satisfy the condition 
that either 0 or 0 Z is prescribed in all regions on the z = 0 plane 
and that 0 is identically zero upstream of a given point. Let 0^ 
and 02 denote two potential functions which satisfy the same boundary 
conditions, and let 0 q denote the potential function formed by taking 
the difference between <f-^ and 02* Mathematically, the potential 
function 0 q is given by 

0 Q (x,y,z) = 0 x (x, y,z) - 0 2 (x,y,z) (53) 

Since 0^ and 0 2 have the same values in certain regions in the 
z = 0 plane then 0q is zero in these regions. Similarly, since 
602 /^z and 60^/ c)z have the same values in the remaining regions of 
the z = 0 plane, then 60 o /oz is zero in these remaining regions. 

The potential function 0q has the boundary conditions that either 0q 
or 60q/6z is zero in all regions of the z = 0 plane and that is 

identically zero upstream of a given point. 

Consider the case where all the boundaries between the regions are 
supersonic. (The slope of the boundaries are such that the component 
of the free stream perpendicular to the boundary is always supersonic.) 

The potential function can be evaluated by use of equations (50) 

and (51) for points in areas which are far enough upstream to be affected 
only by a region where 0 or (f z is prescribed. For all points in 
these areas 0 q is zero as indicated by equations (50) or (51)* It 
follows from equations (50) and (51) that is also zero inside the 

volume above the z = 0 plane, which is affected by these areas alone. 
Thus, the volume where ts identically zero has been moved downstream. 

The same argument can be repeated until the complete z = 0 plane has 
been covered. The preceding arguments cannot be applied to cases where 
the regions have subsonic boundaries) however, if it is permissible to 
distort the boundaries within a strip of infinitesimal width these sub- 
sonic edges can be converted into supersonic edges by replacing every 
element of the subsonic boundaries by a broken line made up of supersonic 
segments. Such a procedure is illustrated in figure 5* If the 
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assumption is made that the subsonic boundaries may be distorted an 
infinitesimal amount, then 0 O is zero over the z = 0 plane and also 
in the space above the z = 0 plane. Equation (53) now reduces to 

0-^(x,y,z) = 0 2 (x,y,z) 

Since 0- L and 0 2 are any two potential functions with the same values 
in certain regions of the z = 0 plane and with the same partial deriva- 
tives with respect to z in the remaining regions, proof has been given 
that only one potential function exists for which the potential is 
prescribed over certain areas in the z = 0 plane and the partial 
derivative with respect to z is prescribed over the remaining areas. 

The boundary conditions for a zero -thickness lifting airfoil with 
a given local angle-of -attack distribution are not of the type discussed 
in the preceding paragraph. The conditions prescribed in the z = 0 plane 
for this type of problem are: The potential function 0 is identically 
zero upstream of the airfoil; 0(x,y,O) is zero except on the plan form 
or in the wake; the partial derivative of the potential with respect 
to z, 0 Z , is given on the plan form; and 0 x (x,y,O) is zero in the 
wake. The preceding boundary conditions do not specify that 0 or 0 Z 
be prescribed in all regions on the z = 0 plane since not 0 but 0 X 
is given in the wake. For airfoils which have trailing edges which are 
always supersonic, the requirement that 0 be continuous in the stream 
direction necessitates the potential in the wake to have the value of 
the potential at the trailing edge of the airfoil. In this case, the 
potential function is uniquely determined. For airfoils which have sub- 
sonic trailing edges the Kutta-Joukowski condition is generally applied 
to the trailing edges to determine 0 uniquely. If the assumption is 
made that the trailing edge can be distorted within a strip of infini- 
tesimal width, then the requirement that 0 be continuous in the stream 
direction can be used to determine 0 uniquely. If the assumption is 
made that the subsonic trailing edge is distorted within the infini- 
tesimal strip so that each segment of each line element of the trailing 
edge is always supersonic (see fig. (6)), then 0 is determined 
uniquely. It is well-known that for airfoils with subsonic trailing 
edges there are an infinite number of solutions which satisfy the 
boundary conditions as stated at the beginning of this paragraph. The 
preceding arguments however prove that there is but one solution for an 
airfoil which has had its subsonic edges replaced by broken lines which 
are always supersonic. Note that it has not been proved that the solu- 
tion obtained by distorting the subsonic trailing edges corresponds to 
the solution satisfying the Kutta-Joukowski condition, nor has it been 
proved that the solution of the distorted trailing edge is independent 
of the manner of distortion. 
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Nonplanar Problems 

The scalar potential resulting from the disturbances caused by a 
nonplanar body can be found from equation (18) provided that both 0 
and S^/dnft are known on some surface Si. Unfortunately, 0 and 
6$/c(nh are not generally known on a surface which fills the requirement 
of the surface Sp; therefore, equation (l8) appears to have little 
value in the calculation of the potential functions for nonplanar systems 
in general. Certain properties of equation (l8) are, however, worth 
investigating. 

The problem of evaluating the potential on the upper surfaces of a 
long rectangular body is discussed. The assumption is made that the 
body extends upstream to infinity and that the sides are parallel to the 
free-stream direction except for small local variations which cause 
small disturbances in the stream. Figure 7(a) shows the forward Mach cone 
from a point on the upper surface of such a rectangular body. This -fig- 
ure also shows that there is a certain part of the surface of the rec- 
tangular body in the forward Mach cone from the point (x,y,z) that can- 
not possibly affect the potential at the point (x,y,z). If the surface 
Si in equation (l8) is taken to be the surface of the rectangular body, 
then equation (l8) indicates that the values of 0 and Sj^dn^ in the 
region which cannot possibly affect the potential at the point (x,y,z) 
should be used in evaluating the potential at the point (x,y,z). The 
only possible explanation of this consideration is that the integral 
of 0 and 80/8n^ caused by the disturbances in the "blind spot” add 
to zero. This consideration can be shown mathematically as follows. 

Let denote the potential function resulting from the disturbances 

inside the blind spot. From equation (U6), it follows that 

\ | 4? s) 45 ' 0 <5>0 

Equation (5^-) indicates that the potential at the point (x,y,z) can be 
evaluated by applying equation (18) to the surface of the rectangular 
body regardless of the blind spots. The same argument holds for other 
bodies with blind spots. 

The preceding arguments can be clarified by a simple illustration 
of the affect of a blind spot. Consider an infinite rectangular body 
such as shown in figure 7 where the only disturbances are caused by a 
small deflected area with a constant slope a with respect to the free- 
stream direction located on the lower surface of the body. The leading 
edge of the deflected area is chosen perpendicular to the free-stream 
direction so that the potential in the region not affected by the 
vertical sides is of a two-dimensional nature. Figure 7(b) illustrates 
such a disturbing surface. 
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by 


The disturbance potential in the two-dimensional region is given 


0(x,z) = 


crv[x - x 2 + p(z - z 2 )J 


where the lower surface of the body lies in the z = z 2 plane and the 
leading edge of the deflected area is in the x = x 2 plane. 


A point on the upper surface which has only the two-dimensional 
flow in its forward Mach cone is illustrated in figure 7(c) . The dis- 
turbance potential for this point is (from equation (l8)) given by 


0(x,y,z) = 2n 


ctV p 2 ffV(| - x 2 )(z - z 2 ) 

~R~ + p d £ dT 1 


Upon performing the indicated integrations the preceding expression 
becomes 

crvjjc - xp - P(z - z 2 )] avfjc - x 2 - p(z - z 2 ) 

23 + 23 

which reduces to 


0(x,y,z) 


0U,y,z) = 0 

This result is a demonstration that the disturbances in blind spots do 
.not contribute to the potential. 

The scalar potential resulting from the disturbance produced by a 
nonplanar body can also be obtained by use of equations (47) and (48) 
provided that the necessary values of 0, 60/6n h , 0', and 80'/8n h 

are known. In dealing with planar bodies 0' could be chosen so that 

0 - 0 ' =0 

and thus equation (47) is reduced to equation (50). Similarly, 0' 
could also be chosen so that 

60 60' 

3% = 5n h* 

and thus equation (47) is reduced to equation (51) • Unfortunately, for 
nonplanar bodies, choosing 0' equal to 0 does not make 60'/6n h * 
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known as was the case for planar problems and, similarly, choosing 
d0 '/ dnj* equal to -80/^% does not make 0’ known. Certain problems 
exist in which 0 ! can be chosen so that 0 can be written as a simple 
integral. 


Intersecting Planes 

Many problems concerning nonplanar bodies deal with disturbances 
produced by two intersecting planes parallel to the free-stream direc- 
tion. In this section, methods of solutions for two planes intersecting 
at various angles are given. The component of velocity normal to the 
surface is assumed to be known. 


Perhaps the simplest case of two intersecting planes occurs when 
the planes intersect at right angles. It is desired to find the poten- 
tial in space resulting from the disturbances produced by the two inter- 
secting planes. This type of problem could represent an isolated cruci- 
form tail with supersonic leading edges undergoing various motions. 
Problems of this type have been solved in references 11 and 12. The 
axes are chosen so that y = 0 and z = 0 are the disturbing planes 
(see fig. 8 ). When y and z are positive, equation (18) becomes 


0 (x,y,z) 



1 afl(s,o + ,p 

R 8 t] 


d 2 f 

L y t 

2 * / 


0(S>o + fS) 


da 


R- 


_ 1 _ 

2ir 



1 d0(S,T),O+) 

R 


da 


P 2 Z t r fl(i,n.o+) 
Js h *3 


da 


(55) 


The surface S]_ has "been taken to he the disturbing surface j thus, S 3 
is the part of the y = 0 plane (z positive) bounded by the z = 6 line 
in the y = 0 plane and the forward Mach cone from the point (x,y,z). 
Similarly, S^ is the part of the z = 0 plane (y positive) bounded by 

the y = 0 line in the z = 0 plane and the trace of the forward Mach 
cone from the point (x,y,z) (see fig. 9 ) • 

The assumption is made that — ^ ^ l!L1 and ) are known 

oil d£ 

and that 0(1, 0 + ,£) and 0(|,q,O + ) are unknown. I'he integrals con- 
taining 0(l,O + ,£) and 0(l,T),O + ) can be eliminated by several appli- 
cations of equation (46). Equation (46) is applied to the volume on 
the left-hand side of the y = 0 plane enclosed by the forward Mach 
cone from the point (x,y,z), the y = 0 plane, the z = 0 plane, and 
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an arbitrary surface upstream of the disturbance (see fig. 10). The 
result of applying equation (1+6) to this volume is 



The surface Sp has been taken to be the y = 0 plane (z positive) and 
the z = 0 plane (y negative); thus, is the part of the 

z = 0 plane in the forward Mach cone from the point (x,y,z) (see fig. 10). 
Adding equations (55) and (56) yields 


00 


r i[30d,o + ,5) 

60'(.^O-,ai 

/ RL di, 

6h 


da - 


P 2 y f 

2 it J 


-^(5,0 + ,C) - 0* ( da - i 


60 ( 1 , 1 - 1 , 0 +) 



R 


da - 


pfz f r 0u,t),q + ) da _ i_ 

R 3 2Tt 


60* (^+1,0+) 


2rt 

pfz f 

2 « , 


R 


da - 


0’(g,O + ,£) 


da 


(57) 


The potential function 0' (£,+!,£) is chosen so that 


0'U,-a,O = 0(l,a,S) 


where a 


is positive. 


In this case. 


a0'(s,o-,o 

611 


60U,Q + ,O 

6ti 
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and equation (57) reduces to 


0(x,y,z) 



1 60(1, 0 + ,£) 

R 6q 


da 


2:r 



lMl2l2£lda 
R 6£ 


P 2 z f r 0(l>n>o + ) 

2 * j s 4 r3 


da 



1 

R 


60 ' (I ,ti,o + ) 
6£ 


da - 



da 


(58) 


Since 0'(l,H,C) is related to 0(1, q,f;) the onl y unknown in the 
preceding equation is 0(i,n,O + ). The region of integration S5 
becomes the part of the z = 0 plane (z negative) in the forward Mach 
cone from the point (x,y,z) obtained by reflecting the disturbing surface 
in the z = 0 plane (z positive) through the z = 0 plane (see 
fig. 10) . 

The problem being considered is one in which the normal derivatives 
of the potential function are known on two planes parallel to the x-axis 
and are intersecting at right angles. The point (x,y,z) has been 
restricted to positive values of y and z. For the present, consider 
the problem of finding the potential above the z = 0 plane when the 
derivative of the potential function with respect to z is known on 
the z = 0 plane, the derivative of the potential function with respect 
to y is discontinuous across the y = 0 plane, and the potential is 
zero a finite distance upstream of a given point. From equation (22), 
the potential for this problem is given by 


0(x,y,z) 




Ml^Mlda 

8t] 



1 60(1, q,0 + ) ^ 
R 6£ 


P 2 z f f 0U,R,O + ) 
2 * J Si r 3 


da 


(59) 


where S 0 is the part of the y = 0 plane above the z = 0 line 
inside the forward Mach cone from the point (x,y,z), and Sq is the 
part of the z = 0 plane inside the forward Mach cone from the 
point (x,y,z). For positive values of y and z equation (59) reduces 
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to equation ( 58 ) if 0 ’ ( 0 + ) is assumed to be the true value 
of 0(t,q,O + ) when q is negative. Since in the original problem 
0(M,9 was not defined when r\ was negative, nothing is violated if 
it is now defined as being l n the region where t] is nega- 

tive and £ is positive. The problem in which the normal derivative 
of the potential function is known on two planes parallel to the x-axis 

and intersecting at right angles has, therefore, been changed to the 

problem in which the derivative of the potential has a known discon- 
tinuity . across the y = 0 plane (z positive) and the normal derivative 
of the potential is known on the z = 0 plane. Note that the poten- 
tial function still remains undefined below the z = 0 plane. 

Since the potential function is undefined below the z = 0 plane, 
it can now be defined so that the resulting potential function is sym- 
metric with respect to the z = 0 plane. Defining the potential below 
the z = 0 plane so that 


0( M,a) = 0( l,n,-a) 

yields the desired symmetry. The result of applying equation (46) to 
the region below the z = 0 plane inside the forward Mach cone from 
the point (x,y,z) is 



where Sg is the part of the y = 0 plane (z negative ) inside the 
forward Mach cone from the point (x,y,z). In applying equation (46) to 
the region below the z = 0 plane, the surface Sq had to be folded 
over the part of the y = 0 plane (z negative) across which 60/6q 
was discontinuous in order to be able to apply equation (46) to this 
region. Equations (58) and (60) can be combined to yield 


0(x,y,z) 


1 [ 1 80(1, 0 + ,g) 1 p 1 80U,q,O + ) 

*,/ s 3 R ^ "*J Sl+ R 8C 



1 S0(|,r 1 ,O + ) 

R 


da - 



1 

R 


80(1, o + ,0 

8q 


da 


1 

TT 


1 

it 


(6i) 
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Equation (6l) contains only integrals of known expressions, and it 
is, therefore, the solution to the problem of two planes intersecting 
at right angles parallel to the x-axis where the normal derivative of 
the potential function is known on both planes. Figure 11 shows the 
cross section of the distribution of velocity normal to the surfaces 
for a problem as represented by equation (55) and its solution as given 


by equation (6l). Note that ^ , 0 __ ) ^ negative) in is the 

reflection of ^ in (t] positive) across the t \ = 0 


plane 


and that 


c>0(l,O + ,Q 


in 


( £ negative) is the reflection of 

in S3 (£ positive) across the £ = 0 plane. This condi- 
tion suggests that the result given by equation (6 1) could also be 
obtained by utilizing the concept of reflecting surfaces. 


The mathematical derivation required for finding solutions to 
problems consisting of two planes parallel to the x-axis intersecting 
at various angles can be reduced by making use of the concept of the 
reflecting surfaces. For this reason, the result given by equation (6l) 
is obtained by use of reflecting surfaces. The potential function can 
be separated into two parts, 0-^ and 0 2 , satisfying the following 
boundary conditions on the disturbing surfaces: 

C>0lU,T),O + ) _ c)0( I ,T) ,0 + ) 

as ' = 

a0 1 (s,o + ,o 

= 0 

d0 2 (l,r),O + ) 
d£ = ° 

= 50(1 0+,g) 

dr) dt] 


A cross section of these boundary conditions is shown in figure 12. 
Only the potential function 0p is treated in detail since the boundary 
conditions for 0^ and 0 2 are of the same type. The normal deriva- 
tive of 0p is zero on the T) = 0 plane; thus, the T) = 0 plane can 
be considered as a reflecting plane. The potential function 0p is, 

therefore, the potential resulting from a distribution of 

which is symmetric with respect to the t) = 0 line and has the value 


,T]>0 + ) 
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¥i(t,T},0 + ) 

of ^ when q is positive. Figure 13 illustrates such a 

distribution. The problem of finding 0j_ has been reduced to a planar 
problem which can be solved by use of equation (50). 

Equation (50) was obtained by defining the potential below the 
z = 0 plane so that the total potential function was symmetric with 
respect to the z = 0 plane. This result caused the derivative of the 
potential function with respect to z to be antisymmetric with respect 
to the z = 0 plane. Figure 14 illustrates the distribution of the 
normal derivative of the function across the z = 0 plane. The problem 
of evaluating the potential function 0p has teen reduced to a planar 
problem. Similarly, the problem of evaluating the potential function 0 2 
can be reduced to a planar problem. Figure 15 illustrates such a pro- 
cedure. The original potential function is the sum of 0-^ and 0g. 
Equation (6l) follows from the preceding results for 0]_ and 0 2 . The 
addition of 0-^ and 0 2 is illustrated by figure l6. 

The concept of reflecting surfaces is now utilized to find the 
potential resulting from two disturbing surfaces parallel to the x-axis 
and intersecting at an angle of 45°. The axis is chosen so that the 
x-axis lies along the intersection of the disturbing surfaces and one 
of the disturbing surfaces lies in the z = 0 plane (see fig. 17). 

The potential function 0 is divided into two parts, 0^ and 0p. The 
boundary conditions on 0-^ and 0 2 are similar to the corresponding 
potential functions used for the disturbing surfaces intersecting at 90°. 
Figure 18 illustrates the boundary conditions for 0^ and 0 2 . The 
surfaces on which the normal derivative of 0]_ is zero can be con- 
sidered as a reflecting surface. This consideration leads to the same 
distribution of the normal derivative of 0^ on the q = 0 plane as 
is given on the £ = 0 plane. Figure 19 illustrates such a distribu- 
tion. The problem of finding 0q for two disturbing surfaces inter- 
secting at A5 0 has been reduced to a problem of two surfaces inter- 
secting at 90°. The solution of 0p can be obtained from equation (6l) . 
Figure 20 shows the surfaces across which the normal derivative of 0^ 
is discontinuous. Since 0]_ and 02 have the same type of boundary 
conditions, then 0 2 has a solution as illustrated in figure 21. The 
original potential function 0 is the sum of 0^ and 0 2 ; therefore, 

0 can be found by considering surfaces of discontinuity as illustrated 
in figure 22. The potential function 0 can be evaluated by use of 
equation (22), because no surfaces across which 0 is discontinuous 
exist and the values of A 60/ckih are known across all surfaces of 
discontinuity. 

Another simple case of two disturbing surfaces parallel to the 
x-axis occurs when the surfaces intersect at an angle of 60°. The 
potential function 0 is divided into two parts, 0^ and 0 2 . The 
boundary conditions on 0^ and 0 2 are similar to the corresponding 
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functions used previously. Figure 23 shows a cross section of these 
boundary conditions. By use of a reflecting surface, the function 0^ 
can be represented by the boundary conditions as shown in figure 24. 

The function 0]_ is undefined for 2k0° of the total angle around the 
x-axis. The function 0p is defined as shown in figure 25. Since no 
surfaces exist across which 0^ is discontinuous, the function 0]_ 
can be evaluated by using equation (22). Similarly, 02 can be defined 
as shown in figure 26. The sum of 0]_ and 02 is illustrated by fig- 
ure 27. The potential function 0 can be found by using equation (22). 

In the preceding paragraphs, methods have been found for deter- 
mining the potential resulting from two plane disturbing surfaces 
parallel to the stream direction intersecting at certain angles. The 
same method can be used to find methods for determining the potential 
resulting from two plane disturbing surfaces intersecting at various 
other angles. 


ROLLING TAILS WITH MULTIPLE RECTANGULAR FINS 


The methods derived in the preceding section are used to find the 
surface velocity potential, the pressure distribution, and the damping 
in roll of rolling tails consisting of four, six, and eight rectangular 
fins. For comparison, these same quantities are also presented for the 
planar tail configurations consisting of one and two rectangular fins. 

An illustration of the tails treated is shown in figure 28. The analysis 
is limited to tail configurations having surfaces of vanishingly small 
thickness and of zero camber. The investigation is also limited to the 
range of Mach numbers for which the region of interference between the 
adjacent fins does not affect the fin tips. 


Tail Consisting of One Fin 

The pressure distribution and the velocity potential on the surface 
of rolling tails made up of one and two rectangular fins can be obtained 
from the results of reference 13. The pressure and potential for the 
tail consisting of only one fin can be found by transforming the axis 
of roll of the tail consisting of two rectangular fins. 

The tail consisting of one fin is divided into regions as shown in 
figure 29(a)* The velocity potential on the surface facing the negative 
y-direction’ is given by: 
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For region I, 
0(x,z) = ^ 

For region II, 


x Z . _! Hi ■ z ) 

T Sln V * 


l/2x 
3(3 " Z 


- b 1/ £ - z 


X«*f-0 


0(x,z) =£f 


(62a) 


(62b) 


For region III, 


0(x,z) = ^ 


+ 0\HR 


(62c) 


For region IV (note that the potential in region IV is the potential in 
region I plus the potential in region III minus the potential in 
region II) , 


0(x,z) = 


XZ . _1 R 

* z . xz 1 M 2 Z ) l/ 2 x \ r/x \ 

p ln v 

X p sin V x + 3 \ P + /V (p ‘ / 




Wz - * + 2^ - 225. 

A 2 p/ 2P 


(62d) 


The pressure -difference coefficient is given by: 
For region I, 


AC = §£_ 

P jrpv 


z sin ' 1 


- 1 + 1 ) 


(63a) 


For region II, 


AC = 

AC p pv 


(63b) 
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For region III, 


For region IV, 


8 p 

" npV 


z sm 


-1 




AC 


P 



z sin - -'- 



f { f - z ) ‘/(I ' # 



( 63 c) 


( 63 d) 


Tail Consisting of Two Fins 

The tail consisting of two fins (fig. 29(b)) has the same potential 
and pressure distribution as a rectangular rolling wing and can, there- 
fore, be obtained from the results of reference 13 . For each tail con- 
sisting of two fins divided into regions as shown in figure 29(b) the 
pressure and potential in regions I and II are the same as the pressure 
and potentials in the corresponding regions for tails consisting of one 
fin. 


Tail Consisting of Four Fins 

Each fin of the tail consisting of four fins is divided into 
regions as shown in figure 29(c). The pressure and potentials in 
regions I and II are the same as the pressure and potentials in the 
corresponding regions for tails consisting of either one or two fins. 

The regions III and IV are affected by the interaction between adjacent 
fins. The potential in region IV is made up of a combination of the 
potentials of regions I, II, and III. Thus, the only real problem is 
the determination of the potential in region III. 

The potential in region III is not affected by the tip and is, 
therefore, the same potential as would be obtained if the fins were 
infinitely long. With the coordinate axes chosen as shown in figure 28(c) 
the point (x,y,z) is restricted to values of y which are negative 
while the values of z are restricted to positive values. Note that 
for a tail with finitely long fins, the potential at a point (x,y,z) 
in the region of interaction is independent of the disturbances produced 
at points located so that their projection on the yz -plane does not lie 
in the second quadrant. The general method previously derived for 
finding the potential resulting from two plane surfaces intersecting at 
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right angles can thus be used to find the potential in the part of the 
region of interaction which is not affected by the tip. 

The velocity component normal to the fin in the z = 0 plane is 
given by 


and the velocity component normal to the fin in the y = 0 plane is 
given by 


Figure 30 illustrates this type of normal -velocity distribution, and 
figure 31 illustrates a cross section of the surfaces of normal -velocity 
discontinuity, which previous results show can be used to obtain the 
potential in the part of the region of interaction which is not 
affected by the tips. 

Note that in figure 32 the discontinuity in the normal velocity 
across the z = 0 plane is the same type as the discontinuity in the 
y = 0 plane. Thus, if the potential resulting from this type of dis- 
continuity (see fig. 32) is known, then the potential resulting from 
any combination of discontinuities of this type can be found. The 
potential for this type of distribution is denoted by 0 O . By the use 
of cylindrical coordinates, as shown in figure 32, the potential at the 
point (x,p,0) can be expressed in terms of 0 q by 


Equation (64) follows from figures 31 and 32. 

The potential function 0 O was evaluated by use of equation (50) 
and is given by: 

For 0 <. 6 <-, 

n; 


0 z (x,y,° + ) = -py 


0y(x,CT,z) = pz 


0(x,p,e) = 0 o (x,p,e) - 0 o (x,p,e - |) 


(64) 


x \l x ^ - P^P 2 + P 2 P^(l - 2 cos2e)in 



+ 




(65a) 
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For 2 — 9 — 


0 Q (x,p,e) = - — 




+ p2p2^ _ g cos 2 0)ln 
pp cos 0 \ 


PP 


X + 


-1 


2 xPp cos 0 sin 

^ ^x 2 - P 2 p 2 sin 2 0 J 
23 2 p 2 sin 0 cos 6 tan~~^ / — — c — — ^ — \ 

\y* 2 - eVj 


(65b) 


From equations (64) and (65), the potential function in the region 
of interaction, which is not affected by the tips, is given by the 

following equation for ^ <6 < n: 


0(x,p,0) = - -2— <23 2 p 2 (l - 2 cos 2 0)ln 


Pp 


\f A 


x + Wx 2 - P 2 p 2 


2xPp [Tos 0 sin 


-1 


pp cos 0 


y\/x 2 - P 2 p 2 sin 2 0/ 


sin 0 sin" 


pp sin 0 


Vx 2 - P 2 p 2 cos 2 0 ; 


oq 2„2 q • 0/4. -1 xlcot 0| _i x|tan 0 

2P p cos 0 sin 0/ tan — ; — - tan — 




|/x 2 - P 2 p 2 


( 66 ) 




The potential in region III of the fin surface is a special case 
= f) P recedin S equation. Thus, the potential in region III 

is given (in Cartesian coordinates) by 
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From equation (67) the pres sure -difference coefficient is found 

to be 


AC p = 


8pz 

xtVP 


sm 


-1 pz 


( 68 ) 


As previously stated, the potential in region IV (see fig. 29) is 
a combination of the potentials in regions I, II, and III. Assume that 
the fins are infinitely long. In this case only two regions (II and III) 
exist since regions I and IV are affected by the tip. The effect of the 
tip can be taken into account by adding a potential which has zero 
normal velocity on the fin and the negative of the pressure of the 
infinite fin in the plane of the fin outboard of the tip. The value of 
such a potential on the fin is given by the difference between the 
potential of region I and the potential of region II. (This potential 
is only the effect of the tip on a semi-infinite rolling wing) . Thus, 
the potential in region TV is the potential in region III plus the dif- 
ference between the potential of region I and the potential of region II. 
Mathematically, the potential in region IV is given by 


0(x,CT,z) = 


xz _i 

T cos 


Pz xz . _1 /p( 2 z ) 

x P I x 

f(f - z - b )\/(I - # + 1 - 1) - z£ 


Pz 


x + 


v'x 2 


P 2 z 2 


(69) 


From equation (69) the pressure -difference coefficient is given by 



Tail Consisting of Six Fins 

The pressure and potentials on the surface of the tail consisting 
of six fins can be obtained in a manner similar to that used for the 
tail consisting of four fins. The pressure and potentials in regions I 
and II are the same as the pressure and potentials in the corresponding 
regions for tails consisting of one, two, or four fins. Regions III 
and IV are affected by the interaction between adjacent fins. The 
potential in region IV is made up of a combination of the potential in 
regions I, II, and III; therefore, the main problem, as for the case of 
four fins, is the determination of the potential in region III. 
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The potential in region III is the same as the potential for a 
tail consisting of six infinitely long fins. The induced velocities 
normal to two of the planes of the fins are illustrated in figure 33* 

For two plane surfaces parallel to the stream direction and inter- 
secting at an angle of 60°, the potential in region III can be obtained 
by a distribution of discontinuities in velocity as illustrated in fig- 
ure 34. Note that the potential in region III can be made up of a com- 
bination of the potentials from a velocity discontinuity as shown in 
figure 35- The potential from this type of discontinuity is denoted 
by 0 O . use of cylindrical coordinates as shown in figure 35 > the 
potential at the point (x,p,0) can be expressed in terms of 0 q by 

0(x,P,6) = 0 o ^x,p,0 - - 0o ( x > p > 0 ' |) " 0o ( x ' p > 0 + 5 ) (?i) 

Equation (71) follows from figures 34 and 35* 

The potential function 0 q was evaluated by use of equation (50) 
and is given by the following equation for 0^0 <_ it: 

0 o (x,p,0) = - P - P - C p ~ — (x - Pp sin 0) (72) 

From equations (71) and (72) the potential function in the region 
of interaction, which is not affected by the tips, is given by the 

following equation for ^ 

0(x,p,0) = £|-j^3"(l ~ 2 cos 2 0) + 2 cos 0 sin 0] (73) 

The potential in region III of the fin surface is a special case 

(0 = ^ of equation (73)* The potential in this region is given (in 
Cartesian coordinates) by 

0 (x,O",z) = (TM 

From equation (74) the pressure -difference coefficient in region III 
is found to be zero. 

The potential in region IV is a combination of the potentials in 
regions I, II, and III) this can be shown in the same way as the poten- 
tial in region IV of the tail consisting of four fins was shown to be a 
combination of potentials from other regions. Specifically, the 
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potential in region IV for the tail consisting of six fins is the 
potential in region III plus the difference between the potentials of 
region I and of region II. Mathematically, the potential in region IV 
is given by 


- 

OJ 

1 N 

* 

— cos- 1 \l 

H 

Tt 

4 

T V 



I(t - 2 - b ) \/(l - z )( 2 + f - 1) 


From equation (75) the pres sure -difference coefficient is given by 



Tail Consisting of Eight Fins 

The pressure and potential on the surface of the tail consisting 
of eight fins can be found by utilizing the potential functions 
used in finding the pressure and potentials on the surface of the tail 
consisting of four and six fins. The pressure and potentials in 
regions I and II are the same as the pressure and potentials for the 
corresponding regions of the other tails. The potentials in regions III, 
IV, V, and VI (see fig. 29(e)) are affected by the interaction between 
adjacent fins. Since the potentials in regions V and VI are combinations 
of the potentials in the remaining regions, the main problem is to find 
the potentials in regions III and IV. 

The potentials in regions III and TV are the same as the potential 
for a tail of eight infinitely long fins. The induced velocity normal 
to two of the planes of the fins is illustrated in figure 36 . From the 
results for two plane surfaces intersecting at an angle of 45 °, the 
potentials in regions III and IV can be obtained by a distribution of 
discontinuities in velocity as illustrated in figure 37. The potential 
resulting from the distribution of discontinuities in velocity as 
illustrated in figure 37 can be obtained from a distribution of dis- 
continuities as possessed by the potential function 0 O used in con- 
nection with the four-finned tail. The potential function used in 

connection with the four -finned tail was evaluated only in the region 
affected by the root sections of the fins. For the case of the eight- 
finned tail, interaction occurs between adjacent fins in regions which 
are not affected by the root sections of the fins. The potential func- 
tion 0 O of the four-finned tail in the region not affected by the root 
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section must be known. In this region, the potential functions 0 q 
for the four- and six-finned tails are the same. 

From figures 32 and 37 the potential function in the region of fin 
interaction which is affected by the root sections of the fins is, 

for g 

0(x,p,e) = -0 o (x,p,e) + $q(x , p ,Q - - 0 O ( X ' P > 0 - I) + tfoif’P ’ 6 + J) 

(77) 

where 0 O is given by equations (65). 

The potential in the part of the region of interaction which is 
affected by the root section is (from equations (65) and (77)), 

for \ 40 ip 


,, „\ 2pp I „ . _i Pp cos 0 

0(x,p ,0) =--^-<-x cos 9 sm 
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^x 2 - P 2 p 2 cos 2 0 
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^x 2 - P 2 p 2 sin 2 (0 + 


Pp sin 0 cos 0 


_1 x cot 0 | . -1 
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^-(l - 2 cos 2 0) 
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x |cot(0 + j)| 
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- + tan" 1 — ! i 141 


Vx 2 - P 2 P 2 


X 2 - P 2 P 2 


(78) 
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The potential in region III is a special case ^0 = of equa- 
tion (78). Setting 0 = ^ in equation (78) yields (in Cartesian 
coordinates) 


0(x,O‘,z) = sin' 1 - x\f2~ sin -1 -j 


pz 


'2x 2 - p 2 z 2 


+ pz tan" 




P 2 z 2 


(79) 


From equation (79) the pressure -difference coefficient is found to 
be given by 


AC p = 


%/sin- 1 »£ - JTsin 

TtSV X 1 


-1 Pz 


JtpV 



(80) 


By inspection of figures (35) and (37) the potential function in 
the region of fin interaction, which is not affected by the root sections 


of the fins, is, for ^ 0 ^ 

0(x,p,e) = 0 o (x,p,0) + 0 o |x,p,0 - 



( 81 ) 

where 0 q is given by equation (72). Substituting equation (72) into 
equation (8l) yields (remember that 0 q is zero upstream of the Mach 

cone from the y-axis), for ^ 0 


0(x,p,e) = - ^)sin 8 + + 

2 

pp 2 cos 0 sin 0 - (l - 2 sin 2 0) - 


^=(cos 0 + sin 0) 
PV2 


Oj if x < ^£-(sin 0 - cos 0) 


V2 


x - ^=(sin 0 - cos 0)j if x > ^£(sin 0 

V2 


^ cos 0 


_ y/2 

Oj if x < Pp sin 0 
x - Pp sin 0; if x > Pp sin 0 


cos 0) 
(82) 
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The potential in region IV is a special case [d = of equa- 
tion (82). Setting 6 = ^ in equation (82) yields (in Cartesian 
coordinates) 


0(x,O",z) = y-[(l - f2)x + 3z] 


(83) 


From equation (83) the pres sure -difference coefficient is found to 
"be given by 


^ ^ 


m 


The potential in region V is the potential in region IV plus the 
difference between the potential in region I and the potential in 
region II j thus, from equations (62a), (62b), and (83), the potential 
in region V is found to be given by 


0(x,O",z) = |jr 




p / 


. .1 Hi - A 

T Bin V x — 


l 2x z b\ 

lf h - z\ 

(z + x 


\T ) \ 

/U ) 

r p 
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(85) 


From equation (85) the pressure -difference coefficient is found to 
be given by 


8p 

= ‘ i(3V 


TtZ 

vf 


- z sm 


+ - z )( z - 1 + 1 ) 


( 86 ) 


The potential in region VI is the potential in region III plus the 
difference between the potential in region I and the potential in 
region IIj thus, from equations (62a), (62b), and (79) > the potential 
in region VI is found to be given by 


0(x,°-,z) = 


-x cos 


-1 pz 


44. - xJ2~ 

x y 


sin 


-1 


(3z 
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( 87 ) 
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From equation (87), the pressure-difference coefficient in 
region VI is found to be given by 



Discussion of Results for Rolling Tails 


Illustrative plots of the chordwise and spanwise pressure distribu- 
tions across one fin for tails with various numbers of fins are shown 
in figure 38* Figure 39 shows illustrative plots of the spanwise loadings 
on one fin for tails made up of various numbers of fins. 

The potential function 0 q used in finding the pressures and 
potentials for the tails consisting of four, six, and eight fins could 
be used in finding pressures and potentials for tails consisting of any 
even number of fins provided that the region of interaction between 
adjacent fins does not affect the tip. The restriction on the region 
of interaction causes the range of validity to decrease as the number 
of fins is increased. The range of validity could be extended, however, 
by use of a pressure or potential cancellation method such as given in 
references 14 and 15 . 

From the potential, the damping in roll per fin was calculated. 

Table I presents the results of these calculations. Figure 40 presents 
the variation of the damping in roll per fin with A3 for tails made 
up of various numbers of fins. For a given Mach number (3 constant), 
figure 40 shows the variation of the damping with aspect ratio. Fig- 
ure 4l presents the variation of the damping in roll per fin with Mach 
number for tails consisting of various numbers of fins with a fin aspect 
ratio of 1.5. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., October 25, 1951 


46 


NACA TN 2641 


REFERENCES 


1. Ursell, F., and Ward, G. N.: On Some General Theorems in the 

Linearized Theory of Compressible Flow. Quarterly Jour. Mech. and 
Appl. Math., vol. Ill, pt. 3, Sept. 1950, pp. 326-348. 

2. Robinson, A.: On Source and Vortex Distributions in the Linearized 

Theory of Steady Supersonic Flow. Rep. No. 9 , College of Aero., 
Cranfield (British), Oct. 1947. 

3. Phillips, H. B. : Vector Analysis. John Wiley & Sons, Inc., 1933* 

4. Hadamard, Jacques: Lectures on Cauchy's Problem in Linear Partial 

Differential Equations. Yale Univ. Press (New Haven), 1923- 

5. Puckett, Allen E. : Supersonic Wave Drag of Thin Airfoils. Jour. 

Aero. Sci., vol. 13, no. 9, Sept. 1946, pp. 475-484. 

6. Heaslet, Max. A., and Lomax, Harvard: The Use of Source-Sink and 

Doublet Distributions Extended to the Solution of Boundary -Value 
Problems in Supersonic Flow. NACA Rep. 900, 1948. (Formerly 
NACA TN 1515.) 

7. Martin, J. C.: Retarded Potentials of Supersonic Flow. Quarterly 

Applied Math., vol. VIII, no. 4, Jan. 1951, PP* 358-364. 

8. Kellogg, Oliver Dimon: Foundations of Potential Theory. Frederick 

Ungar Publishing Co. (New York), 1929- 

9. Miller, Frederic H. : Partial Differential Equations. John Wiley & 

Sons, Inc., 1941. 

10. Lomax, Harvard, Heaslet, Max. A., and Fuller, Franklyn B. : Formulas 

for Source, Doublet, and Vortex Distributions in Supersonic Wing 
Theory. NACA TN 2252, 1950* 

11. Bleviss, Zegmund 0.: Some Roll Characteristics of Cruciform Delta 

Wings at Supersonic Speeds. Preprint of paper presented at 

Annual Summer Meeting, Inst. Aero. Sci. (Los Angeles), July 12, 1950. 

12. Ribner, Herbert S.: Damping in Roll of Cruciform and Some Related 

Delta Wings at Supersonic Speeds. NACA TN 2285, 1951. 

13. Harmon, Sidney M.: Stability Derivatives at Supersonic Speeds of 

Thin Rectangular Wings with Diagonals ahead of Tip Mach Lines. 

NACA Rep. 925, 1949* (Formerly NACA TO 1706.) 


NACA TN 2641 


47 


14. Goodman, Theodore R. : The Lift Distribution on Conical and Non- 

conical Flow Regions of Thin Finite Wings in a Supersonic Stream. 
Jour. Aero. Sci., vol. 16, no. 6, June 1949, PP- 365-374. 

15. Mirels, Harold: Lift-Cancellation Technique in Linearized 

Supersonic -Wing Theory. NACA TN 2145, 1950. 


48 


NACA TN 2641 


TABLE I 

DAMPING- IN-ROLL COEFFICIENT PER FIN 


Number of 
fins 

j3Ci p per fin 

Valid. 

for 

1 

1 + 4A P - 24A 2 p 2 + 32A 3 p 3 
24 a 3 P 3 

A P ;> r 

2 

l + 8A0 - 48A 2 p 2 + 64 a 3 P 3 
48A 3 p 3 

H|C\J 

All 

ca 

< 

4 

4 

*A3p3 

-i + -g^(i + 8A P - 48 a 2 P 2 + 64a 3 P 3^ 
9 19*- \ j 


A(B > 1 

6 

4 

A3p3 

- — . + — - (1 + 8A P - 48a 2 (3 2 + 64a3 P 3J 
9^3 192 V L 


»*is 
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_ ^ [1 ^ - + -^-(l + 8 a P - 48A 2 (B 2 + 64a3 P 3) 

a 3 3 3 L 9* 192 V 'J 

A P > (2 
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figure 1.— Cross section of the region of integration used in connection 

with equation (9). 
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Figure 2.- Cross section of the region of integration used in connection 

with equation (lit) • 
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Figure 3. A cross section of the region in the forward Mach cone from 
the point (x,y,z) showing the surface Si enveloping a surface of 
discontinuity S. 
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Figure 4.- Cross section of the region of integration used in connection 

with the potential function 0' . 
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(a) Subsonic element within 
infinitesimal strip. 


(b) Broken line replacing 
line element. 


Figure 5- _ An illustration of a method of replacing a subsonic line ele- 
ment with a broken line made up of supersonic segments. 



(a) Wing with subsonic edge. (b) Line element at a break in a 

subsonic trailing edge. 



(c) Possible distortion of a subsonic trailing edge at a break in a sub- 
sonic trailing edge. 

Figure 6.- Methods of distorting a subsonic trailing edge to determine the 

potential function uniquely. 
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Forward 
Mach cone 


Intersection of the 
forward Mach cone 
with the rectangular 
body 


Z 




(a) The forward Mach cone from a point on the upper surface of a rectangular 

body. 


Figure 7* — Rectangular body parallel to free— stream direction. 
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(b) Bottom view of rectangular body with deflected area on the lower 

surface. 


FORWARD MACH CONE 



(c) Side view of rectangular body and the forward Mach cone from the 

point (x,y,z). 


Figure 7 .- Concluded. 
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Figure 9 .- Regions of integration for equation (55). 
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Figure 10.- Regions of integration for equation (56). 
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(b) Solution to problem as given by equation (6l). 

11.— A cross section of the distribution of the velocity component 

normal to the z = 0 and the y = 0 planes represented by equations (65) 
and (61). 


6o 


NACA TN 2641 


0 — 0\+ 02 
TTm-rTTfit i trrTTrrrn. 


01 

nmmnnmn. 



Figure 12.- A cross section of the distribution of velocity normal to 
the disturbing surfaces for the potential functions 0, 0-^, and 02- 
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Figure 13. - The reduction of 0]_ to a planar problem. 
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Figure 14. - The normal derivative of the potential function across the 
z = 0 plane obtained by applying equation (50) to a planar problem. 
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Figure • An illustration of the reduction of 02 to a planar problem. 
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Figure 16.- An illustration of the addition of 0^ and 02* 
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Figure 17.- Position of coordinate axes for disturbing surfaces 

intersecting at 1*5°. 





.figure 18. A cross section of the velocity distribution on the disturbing 
surfaces for the functions 0]^ and 02 for the disturbing surfaces 
intersecting at kS° . 
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Figure 19.- Reflection of the normal derivative of 0p on the 

y = 0 plane. 




Figure 20.- A cross section of the surfaces of discontinuity which can 
be used to evaluate the potential function 0^ f° r two surfaces 
intersecting at 45°. 
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Figure 21.- A cross section of the surfaces of discontinuity which can 
be used to evaluate the potential function 0 2 for two planes 
intersecting at U5°. 



Figure 22.- The addition of 0 X and 0 2 to obtain the potential func- 
tion 0 for two surfaces intersecting at h5° . 
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Figure 23.- A cross section of the boundary conditions of the functions 0-^ 
and 02 for two disturbing surfaces intersecting at 60°. 




Figure 2k. - The reflection of 0-j. through one surface. 
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Figure 25.- A method of defining so as to eliminate discontinuities 

in the potential function. 




Figure 26.- A method of defining so as to eliminate discontinuities 

in the potential function. 
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Figure 28.- Types of tails treated. 
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(a) One fin. 


(b) Two fins, 



(c) Four fins. 
Z 




\ 1 


/ \ 

/ HL 


adjacent f/n 


(d) Six fins. 



/ Mach /me from 
A — adjacent f/n 

/ March /me from 
root sect/ on 


X 


^X 


Mach lines 


(e) Eight fins. 

Figure 29.- Regions of similar disturbances for tails consisting of 

rectangular fins. 
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Figure 30.- A cross section of the 
normal-velocity distribution on 
two plane surfaces representing 
the region of interaction for a 
rolling tail with four fins. 



Figure 31 •- A cross section of the 
velocity discontinuity distribu- 
tion used to find the potential 
in part of the region of inter- 
action for a tail of four fins. 



Figure 32.- A cross section of the velocity discontinuity distribution 
associated with the function 0 q for a tail with four fins. 
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Figure 33. - -4 cross section of the 
velocity distribution normal to 
the planes of two fins of a 
rolling tail consisting of six 
fins. 



Figure 34.- A cross section of the 
velocity discontinuity distribu- 
tion used to find the potential 
in part of the region of inter- 
action for a tail of six fins. 



Figure 35>.- A cross section of the velocity discontinuity distribution 
associated with the function 0 q for a tail of six fins. 
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Figure 36 .- A cross section of the normal velocity induced on the planes 
of two fins of a rolling tail consisting of eight fins. 



Figure 37 A cross section of the velocity discontinuities which can be 
used to obtain the potential in the region of interaction between 
adjacent fins for a tail consisting of eight fins. 
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(a) Tail of one fin. 

Figure 38 .- Chordwise and spanwise pressure distributions on a fin of 
aspect ratio 1.5 at a Mach number of for tails consisting of 

rectangular fins. 
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(b) Tail of two fins. 
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Figure 38 .- Continued 
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(c) Tail of four fins 
Figure 38 .- Continued 
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(d) Tail of six fins. 
Figure 38 .- Continued. 
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(e) Tail of eight fins 
Figure 38.- Concluded. 
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Figure 39.“ Spanvri.se loading on a fin of aspect ratio 1.5 at a Mach number 
of \fZ for tails consisting of rectangular fins. 
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Figure 40.- The damping in roll per fin for tails consisting of rectangular 

fins. 
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Figure Ul. An illustrative variation of per fin with Mach number 

for tails consisting of rectangular fins with a fin aspect ratio of 1.5. 
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